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Abstract 

Some of the basic concepts regarding asymptotic series are re- 
viewed. A heuristic proof is given that the divergent QCD pertur- 
bation series is asymptotic. By treating it as an asymptotic expansion 
we show that it makes sense to keep only the first few terms. The ex- 
ample of e^e~ annihilation is considered. It is shown that by keeping 
only the first few terms one can get within a per cent (or smaller) of 
the complete sum of the series even at very low momenta where the 
coupling is large. More generally, this affords an explanation of the 
phenomena of precocious scaling and why keeping only leading order 
corrections generally works so well. 
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By virtue of its property of asymptotic freedom QCD perturbation the- 
ory gives an excellent description of many high energy processes involving 
large momentum transfers, especially in inclusive phenomena. The agree- 
ment with theory often remains valid down to surprisingly low momenta as 
in the classic case of precocious scaling in deep inelastic lepton scattering 
where approximate scaling persists down to less than a GeV. At first sight 
the great success of perturbation theory is all the more surprising given its 
relatively large coupling constant (g) and the fact that the series is divergent. 
Furthermore, it is generally believed that the nature of the divergence is suf- 
ficiently severe that the series is not summable by conventional methods: no 
available technique exists for reconstructing a unique analytic representation 
for its sum. This is in marked contrast to scalar field theories which have 
been shown to be Borel summable and for which a well-defined representa- 
tion exists [|T|. Thus a major theoretical question concerning QCD is how to 
control and make sense out of the non-summable divergence in its perturba- 
tion expansion. Ultimately large order loop contributions must dominate the 
leading terms so a well-defined procedure justifying the use of just the leading 
order estimates is required. Equally as important is to understand the errors 
incurred by such a procedure. A natural framework for dealing with this 
problem is to treat the QCD perturbation series as asymptotic and apply an 
Euler-Poincare analysis to estimate the number of terms that should be kept 
in order to obtain an optimal estimate of its sump|. A by-product of such an 
analysis is an estimate of the error. In what follows we first give a heuristic 
proof that the series is indeed asymptotic. Using generic forms for estimates 
of the large order coefficients optimal sums for the series are subsequently 
obtained. It is worth emphasizing that the presumption of an asymptotic 
series is significantly weaker than assumptions required for summability. 

Generally speaking, the nature of an asymptotic series is such that, as the 
first terms are calculated, the correct sum is uniformly approached until, 
after an optimum number of terms, Nq, is reached, adding additional terms 
drives the sum further and further from the correct result ultimately leading 
to a divergence. A crucial characteristic of this behaviour is that by retaining 
only the optimal A'o terms and discarding the rest, one can get exponentially 
close to the correct sum when the expansion parameter is small. We shall 
show that carrying through such an analysis for QCD leads to the following 
conclusions: 
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i) No is a relatively small number which depends on the characteristic 
momentum. 

ii) At high energies, where as{= g^/An) ~ 0.1, the error incurred by 
keeping only Nq terms is much less than a per cent. 

iii) At lower energies, where as ~ 0.2, say, Nq decreases and can become 
less than 2; the corresponding error can still be less than a per cent. 

Thus, at low energies where as is quite large, the sum of the perturbative 
series can be approximated to within a per cent by keeping only the first 
couple of terms! Adding higher order loop contributions will only drive one 
further from the correct result. Since the error is only roughly a per cent or 
so, it is still well within typical experimental errors. This, therefore, offers 
a possible explanation as to why perturbation theory works so well even at 
rather modest energies and momenta where as is not so small. The nature of 
these conclusions is rather general (when applied to appropriate processes); 
the details, however, will be process dependent and, in general, depend on 
the effective expansion parameter. 

Before deriving these results we first review the definition and some salient 
properties of asymptotic series. Consider a function F{z) which is analytic 
in a wedge centered on the origin. If the wedge opening is less than 27r, then 
F{z) is non-analytic at the origin and so a power series development 

oo 

F{z)^Y.^nZ'' (1) 


must have zero radius of convergence. For QCD perturbation theory, z = g"^ 
and the the sum of all Feynman graphs at order n. For fixed z the 

remainder 

N 

RN{z) = \F{z)-Y.a„z^\ (2) 



diverges for large N. On the other hand, if the series is asymptotic, then, for 
a given A^, Rn{z)/z^ ~ when 2; ~ 0. This follows from the basic definition 
of an asymptotic series[0], namely that, for a fixed A^ and z ~ 0, 

Rn{z) < C^z""^' (3) 
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Now let us exploit, the analytic structure of F{z) by deriving a general- 
ized dispersion relation. Consider some point z encircled by a contour C 
in a region where F{z) is everywhere analytic and write a standard Cauchy 
representation 

_ f (4) 

^ ^ Jc 2mz' - z ^ ^ 

The contour can be distorted to encircle all the singularities of F[z) leading 
to 

F{z) = f (5) 

where f{z) = {2txi)^^ X^disc F{z) with the sum and integral being taken over 
all discontnuities of F. Possible contributions from the contour at infinity 
have been dropped; if the integrals are not sufficiently convergent for this 
to be valid, sufficient subtractions are assumed to have been made. Their 
presence does not change the general argument. Formally expanding (|^) in 
powers of z leads to 

m (6) 



L 2"+^ 



This is the basic formula used to derive the a„. Indeed, if a path integral 
representation is used for F{z) then this formula reproduces conventional 
Feynman graphs. It can also be used to estimate the large n behavior of a„ 
via a saddle point technique. The path integral representation for F{z) has 
the generic structure / D^e"'^''^'/^, where S is the action; [here the generic 
field stands for all fields including fermions and gauge bosons]. Thus, 
singularities in z develop only when Re z < 0, i.e., along the negative real axis 
in which case f{z) reduces to 2i lmF{z)^. This relatively simple structure 
is broken by renormalization which gives rise to a considerably more complex 
singularity structure This situation makes it difficult to give a rigorous 
proof that the series must be asymptotic; nevertheless, a heuristic proof can 
be given. 

The crucial observation is that Eqs. @ and (|^) can be combined to give 
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If, when z ~ 0, the integral exists, then the constraint of Eq.(^) is clearly 
satisfied; however, in this limit the integral is, at least formally, simply Qn+i 
[see Eq. (^]. Furthermore, these coefficients, a„, must exist, since they 
represent the sum of all Feynman graphs of a given order, and can, in fact, 
be generated from Eq. (||). Thus, as 2; ~ 

Rr,{z) ~ a;v+/^^+' (8) 

which manifestly satisfies thereby showing that the series must be asymp- 
totic §. Q.E.D. 

The expression manifests one of the essential features of asymptotic 
series: as increases, grows whereas diminishes, thereby leading to a 
minimum value Rn- Before applying this to QCD let us first briefiy review 
the relationship of these ideas to summability. 

The idea behind Borel summability is that, although the series (|I]) may 
be divergent, the series 

G{z)^t^-^z- (9) 
"'■ 

may be convergent. In that case F{z) can formally be obtained from G(z) 
by the Laplace transform 

Fiz) = - / rfMe^"/^G(M) (10) 
z Jo 

This then serves to define a function that is asymptotic to the original diver- 
gent series. Under certain restrictive conditions (loosely speaking, that the 
integral and its transform exist) this construction defines a unique function. 

As already mentioned it is believed that the Borel procedure can be ap- 
plied to scalar field theories to give a representation for its perturbative sum; 
such is not the case, however, for QCD. The difference between the two 
can be characterized by the classic examples (i) a„ = {—l)""n\ for which 
G{z) = (1 + z)~^ and (ii) a„ = n\ for which G{z) = (1 — z)^^. The latter, 
which is the analog of QCD, is not Borel summmable since the integral is 
not well-defined because of the ambiguity in how to treat the singularity at 
M = iP|. The former, however, which is the analog to scalar field theory, is 
well-defined. From Eq. (^O]) we deduce that 

E{z)=r^ (11) 

^ ' Jo l + vz ^ ^ 
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is the analytic reconstruction of the divergent series 

I- z + 2\z'^ -?,\z^ + A\z^ - (12) 

Eq. (p!T| ) can therefore be used to define and evaluate the divergent sum, Eq. 

(|12D. Thus, for example, a numerical evaluation gives £'(0.1) = 0.915633 

However, and this is the important point, it is also possible and, in general, 
considerably easier and more efficient to use the series directly to get an 
excellent estimate for E{z). To see how this comes about, we return to 
Eq. (H) from which we see that, for small z, Rn{z) ~ (A^ + l)!z^^^. Thus 
Rn first decreases as increases, reaching a minimum at = A^o? where 
dRN{z) I dN = 0, and then increases in an unbounded fashion. The minimum 
occurs at A^ = A^o ~ {l/z)e^^'^^ — 2 at which point 

Rno ~ [27ie{No + 2)]i/2g-(iVo+2) ^ (2vr/^) VSg-V- (13) 

This is a remarkable result, for it shows that one can get exponentially close 
to the correct answer by keeping only the first A^o terms. Eq. ([T3| ) thus 
represents the closest one can get to the correct answer by this technique. So, 
if z = 0.1, as in the above example, this gives A^o ~ 8.5. Keeping 8 terms in 
(|12|) then gives E(O.l) ^ 0.915819 whilst 9 terms gives E{0.1) ^ 0.915460 in 
excellent agreement with the exact result obtained by numerically evaluating 
the integral. The error thus incurred is in agreement with Eq. (pTsP which 
gives i?7Vo ~ 3.6 x 10~^. Now, suppose that z = 0.2, then A'q ~ 5.5 so only 
3-4 terms need be kept! The error thus generated is thus only R]\f^{0.2) ^ 
3.78 X 10^^. Keeping 3 or 4 terms one trivially finds that, E{0.2) lies between 
0.832 and 0.870 to be compared with the exact, number 0.8521 ... Thus, if 
a given accuracy is sufficient one need only keep a relatively small number 
of terms in the series. Even when z ^ 0.2 one can get within a few per cent 
of the exact answer. In such a situation, there is no advantage in having 
an explicit representation such as (p^ and little or nothing would be lost 
without it. We now apply these ideas to QCD and show how the asymptotic 
nature of the series and the large value of can be put to advantage. 

As an illustrative example consider the QCD contribution to the total 
e~^e~ annihilation cross-section 0. This is usually expressed as a ratio, R, 
normalized to the cross-section for fi'^^^ production. Asymptotic freedom 
dictates that its high energy behaviour is ultimately given by J2Qh Qi be- 
ing the charge of the ith quark species. Leading corrections to this can be 
developed in a standard perturbative series: 
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R 



-.2 



oo / ^2 



EQljEri^K^) (14) 



Here q is the four- momentum delivered by the e~^e~ pair; we have also 
made manifest the arbitrary scale /i needed to define as- The first four 
coefficients of this series have been calculated: when = /i^ they are 
1^0 = f'l = ^',f2 = 1-41 and = —12.8. The first two are both scheme 
and momentum independent whereas the numbers quoted for r2 and are 
in the modified-minimal-subtraction [MS) scheme with five flavors. 

It is generally agreed that for large n the structure of the r„ has the 
generic form 

r„(l) ^ CaVr(n + c) (15) 

where C, a, b and c are constants. This can be motivated from the path 
integral representation where all of the coupling constant dependence resides 
in the factor e~^^^ . Note that this implies that the real effective expansion 
parameter is aas/n rather than a^/vr. From Eq. (|15|) we can also estimate 
the optimum number of terms to be kept in order to get as close as possible 
to the exact sum: 



No^a-^ -1-c (16) 
as well as the corresponding error: 

Rn^ ^ (27r)i/2^a^/2-^-'=e-i/'^ (17) 

Now, renormalisation, which introduces the arbitrary scale fi into the 
problem, forces i? to be a function of the single variable z = (g^/ /U^)e~ -f '^^/Pia) 
rather than of the two variables j ^ and as(/i) separately. Here iii^g) is the 
usual /3-function which has the perturbative expansion: i^i^g) ~ —g^{h\ + 
^29^ + ...). Thus z ~ (g^//i^)e^/''i^ g^' where h' = h2/h\. This suggests that 
the real expansion parameter is not as/vr but rather hig^ = Anbias since these 
must always occur together. We have argued [0 from a more detailed analysis 
based on analyticity and the use of a saddle point technique that it is, in 
fact, An'^ebias so that, in terms of Eq. (|15]), a = in^ebi, b = —2 and c = b'. 
The general structure of the results and conclusions presented here do not 
depend on the detailed values of the coefficients. Although there has been 
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some criticism of the derivation of the estimates in ||^ it is the generic 
structure represented by Eq. (|^) with a ~ in^bi that drives the general 
conclusions. Note also that an expansion around instanton solutions, where 
S Stt^ would suggest a value of a = 1/2 which is significantly smaller 
than the value derived from the above renormalisation argument. These 
latter "non-perturbative" effects ( "renormalons" ) can therefore be expected 
to dominate the large n behaviour arising from instantons. 

We can therefore estimate that No ~ (in'^ebias)'^ — b' — 1. With as ~ 

0. 15, this gives No ~ 5. Thus it makes sense to keep only 5 terms at best 
in the series. Keeping more than this, even if calculable, drives one further 
from the correct sum. An estimate of C which gives good agreement with 
r3(l) is (— 47r^6ie~^ With this the miminum error is found to be 7 x 10~^, 

1. e. less than 0.1%. This is illustrated in Table I and Fig.l. 

The discussion so far has focused on the single value = /x^. Indeed, all 
of the numbers quoted are scale dependent, the scale being determined by 
the value of n for which as = 0.15. It is conventional to use the invariance 
of R to changes in to transfer the dependence from the to as by 
introducing a running coupling constant: 



R 



-^,OLs{^i) 



R 



l,as{q^ 



as 



IT 



(18) 



where J^^ das/ P{as) = ln(g^/yU^) with /3{as) = gP{g)/167i^. This is usually 
expressed in terms of a QCD scale parameter A: as(g^) = (47r6it)~^[l — 
b' Int/t + ■ ■ ■], where t = ln(g^/A^). The analysis that was applied to the 
original series in terms of asifJ^) can now be applied to Eq. ([T8|). The only 
difference is that as is replaced by dsiQ"^) so that the corresponding No now 
becomes momentum dependent. From Eq. (|TB|) we obtain an expression for 
the -dependence of the optimum number of terms for the series (KT 



No ^ e-\t + b' Int) + 1/2 - b' 



(19) 



For simplicity we have kept here only the leading term; corrections have 
only a small effect. The corresponding error at this optimum number is 
given approximately by Rno{q'^) ~ {87r^)~^^'^bi^ (A"^ / q'^Y^^ . These equations 
exhibit the behavior already alluded to, namely that, as decreases and 
as(g^) increases, the optimum number of terms that can be retained actually 
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Figure 1: Plot of logio Rn-i vs. for various values of with A = lOOMeV. 
The curves are labeled by both the values of in GeV/c^ and the equivalent 
value of as- Notice that the error in keeping only the one-loop correction 
{N — 2) is always < 1% even at very low values of q^. The locus and 
minimum value are approximately given in the text. 

decreases. The price to be paid for this remarkable behaviour is that the 
error thus incurred correspondingly increases. 

This pattern is shown in Fig. 1 where the error Rn is plotted versus for 
various values of q'^ with A = lOOMeV. As can readily be seen from the graph 
only a very few terms need be retained in order to get within a few per cent of 
the correct sum even down to very low values of g^. Explicit non-perturbative 
corrections to the sum arising, for example, from instantons are expected 
to be of order e"^'^/"" which is much smaller than even the smallest value 
oi Rn. This call be put shghtly differently by noting that such instanton 
contributions behave like (A^/g^)"^'^ which remains much smaller than 
RNo{(f) down to values of q^ ~ IGeV/c^. This therefore serves as a possible 
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explanation as to why tree graphs, supplemented by one loop corrections, give 
such a good description of many processes in QCD even at rather modest 
energies where is relatively large. Both higher order terms and explicit 
non-perturbative contributions typically contribute less than a per cent; only 
when one is sensitive to infrared, bound state or threshold problems is it 
expected that the perturbative description becomes inappropriate. 

Finally we should draw attention to the enormous size of the coefficients 
ttv in Table 1; eventually these lead to contributions that exceed the leading 
term which, if taken seriously, would invalidate the predictions of asymptotic 
freedom as is clear from Fig. 1. The technique suggested here says that, 
in the spirit of asymptotic series, all contributions from diagrams of 0{No) 
and higher should be ignored. It could therefore be argued that, just as in 
the series (|12|) where an individual term such as 100!(0.1)^°'^, for example, 
should be thought of as irrelevant and meaningless as far as contributing to 
the " sum" of the series is concerned, so in ([1^ 8th order graphs, for example, 
are likewise meaningless and irrelevant (at least until one gets to sufficiently 
high energies). 
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Table 1: The coefficients occurring in the expansion of R, Eq. (|14D- The 
ffist four are exact while the remaining are estimates from Eq. (|l^). Also 
shown are estimates of Rn, the corresponding deviation from the exact result 
of the sum of the ffist N terms of the expansion (calculated using as = 0.15). 
Note that although the rapidly diverge, Rn goes through a minimum near 
~ 5 where the error is only 6.85 x 10~^. 



N 


rN+i 


Rr, = \rr,+iiaj7ry'+'\ 


-1 


1 


1 





1 


4.77 X 10-2 


1 


1.41 


3.21 X 10-3 


2 


12.8 


1.39 X 10-3 


3 


163.4 


8.49 X 10-^ 


4 


2.7 X 10^ 


6.71 X 10-^ 


5 


5.78 X 10^ 


6.85 X 10-^ 


6 


1.53 X 10^ 


8.63 X 10-^ 


7 


4.78 X 10^ 


1.29 X 10-3 


8 


1.74 X 10^ 


2.25 X 10-3 


9 


7.23 X 10^° 


4.45 X 10-3 


10 


3.38 X 10^^ 


9.93 X 10-3 


11 


1.75 X 10^^ 


2.45 X 10-2 


12 


9.97 X 10^^ 


6.68 X 10-2 


13 


6.19 X lO^*^ 


1.98 X 10-1 


14 


4.17 X 10^8 


6.38 X 10-1 


15 


3.02 X 10^1 


2.2 
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